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Abstract. Based on the quasi-particle description of the QCD medium at finite temperature and density 
we formulate the phenomenological model for the equation of state that exhibits crossover or the first 
order deconfinement phase transition. The models are constructed in such a way to be thermodynamically 
consistent and to satisfy the properties of the ground state nuclear matter comply with constraints from 
intermediate heavy-ion collision data. Our equations of states show quite reasonable agreement with the 
recent lattice findings on temperature and baryon chemical potential dependence of relevant thermody- 
namical quantities in the parameter range covering both the hadronic and quark-gluon sectors. The model 
predictions on the isentropic trajectories in the phase diagram are shown to be consistent with the recent 
lattice results. Our nuclear equations of states are to be considered as an input to the dynamical models 
describing the production and the time evolution of a thermalized medium created in heavy ion collisions 
in a broad energy range from SIS up to LHC. 

PACS. 21.65.+f Nuclear matter - 24.85.+p Quarks, gluons, and QCD in nuclei and nuclear processes - 
12. 38. Aw General properties of QCD - 12.38.Mh Quark-gluon plasma 
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1 Introduction 

QCD at the finite temperature T and/or baryon chem- 
ical potential fis is of fundamental importance, since it 
describes the relevant features of particle physics in the 
early universe, in neutron stars and in heavy-ion colli- 
sions (see e.g. 0E]). With the relativistic heavy ion col- 
lision experiments at AGS, SPS and RHIC accelerators 
one explores the phase diagram of strongly interacting 
matter in a broad parameter range of temperature and 
baryon density. Lattice QCD results on the Equation of 
State (EoS) of QCD matter provide a basic input for the 
analysis of experimental signatures of a possible quark- 
gluon plasma formation in heavy-ion collisions. Directly 
addressing the EoS, hydrodynamics realizes the connec- 
tion between the matter properties and observables. The 
hydrodynamic treatment of the whole time-evolution of 
colliding nuclei requires knowledge of the nuclear EoS 
within a large interval of its thermodynamic variables cov- 
ering both quark-gluon and hadronic sectors. 

In the recent years a significant progress has been made 
in understanding the phase diagram of QCD at non-zero 
baryon chemical potential as the nonperturbative lattice 
QCD methods were extended to access the relevant re- 
gions of the phase diagram. Recently, the first lattice cal- 
culations have been performed for a non-vanishing T and 
HB for systems with N f = 2 [3j and N f = 2 + 1 fla- 
vors. However, due to a set of approximation the Lattice 



Gauge Theory (LGT) is still not able to provide results 
on the properties of the hadronic matter in the confined 
phase. LGT is also restricted to moderate values of the 
baryon chemical potential (1b such that \ib < T. That is 
why different phenomenological models are required to de- 
scribe thermodynamic properties and equation of state of 
QCD matter for larger baryon densities. Obviously, such 
models depend on the set of parameters that are usu- 
ally fixed to reproduce existing LGT results as well as 
the basic phenomenological properties of the nuclear mat- 
ter obtained from the experimental data. Recently, the 
thermodynamics of the quark-gluon phase was interpreted 
quite successfully within the QCD inspired massive quasi- 
particle models |^|7|$B ?|10I11I12I13I14| . On the other hand, 
demonstrated by lattice calculations, a rapid growth of the 
energy density e and pressure p when approaching the crit- 
ical temperature T c was shown to be reproduced in terms 
of the hadron resonance gas model with scaled masses ^1 
116) . Only recently there were attempts to describe lattice 
QCD thermodynamics both above and below T c in terms 
of a field theoretical model, including features of both de- 
confinement and chiral symmetry restoration |17| . as well 
as within some phenomenological models that are based 
on lattice QCD results for the quark-gluon partition func- 
tion |18| . Some unique parametrization of the QCD EoS 
below and above T c was also presented in Ref. |19j . 

The phenomenological equation of state should be not 
only thermodynamically consistent |2()| but should also 
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be capable to reproduce the global behavior of the nuclear 
matter near the ground state and its saturation properties. 
In addition, there are experimental restrictions coming 
from the flow analysis in heavy-ion collisions which limit 
the acceptable theoretical values of pressure in a finite 
interval of baryon densities tib at T = [S"Q. Some con- 
straints on the EoS are also imposed through the analysis 
of cold charge-neutral baryonic matter in /3-equilibrium 
compact stars |22II23| . There are also essential constraints 
on the model properties coming from the recent LGT re- 
sults. 

In this paper, we will construct the EoS of a strongly 
interacting QCD matter with deconfincmcnt phase tran- 
sition that satisfies the above mentioned hadronic con- 
straints and those imposed by the recent lattice QCD re- 
sults obtained for (2 + 1) - flavor system at the finite T 
and non-vanishing baryonic chemical potential. 

The paper is organized as follows: In Section 2 we in- 
troduce the quasi-particle model for the EoS with decon- 
finement phase transition. In Section 3 the model predic- 
tions are compared with the recent lattice data obtained in 
(2+l)-flavor QCD at the finite T and ps- Our results and 
comments on the properties of the QCD equation of state 
and thermodynamics are summarized in the last Section. 



2 The Equation of State 



Lattice results show that even at temperature T much 
larger than deconfincmcnt temperature T c the thermo- 
dynamical observables like pressure or entropy, baryon 
number and energy density are still by > 20% deviating 
from their asymptotic ideal gas values. Such deviations 
observed at T > 2T C were shown to be well understood 
by a systematic contribution in the self-consistent imple- 
mentation of quasi-particle masses in the HTL-resummcd 
perturbative QCD [23 . On the other hand, the LGT ther- 
modynamics below T c was shown to be well reproduced 
by the hadron resonance gas partition function |15Ulb |. To 
possibly describe the thermodynamics at T — T c or near 
the phase transition an additional model assumptions are 
required [2*5115] . 

It is clear from the above that the straightforward 
model for the QCD EoS can be constructed by connecting 
a non-interacting hadron resonance gas in the low tem- 
perature phase with an ideal quark gluon-plasma in some 
non-perturbative bag for the color deconfined phase [H"] . 
These phases are matched at the phase transition bound- 
ary by means of the Gibbs phase equilibrium condition. 
By construction, this approach yields the first order phase 
transition. Such MIT bag-like model ^26] is so far the sim- 
plest method to implement the confinement phenomenon 
in the EoS, though it has some serious shortcomings. 

A more complete method to model QCD EoS is based 
on the effective Hamiltonian that includes interactions of 
the constituents. In the quasi-particle approximation such 
Hamiltonian can be modelled through density-dependent 



mean-field interactions |20U27U28| : 
H = J2 E/ dr^(r,s) ■ 

jeh,q,g s 

( ^-V2 + mf + U 3 {p) \ ^(r, s) - C{p) ■ V , (1) 

where j enumerates different species of quasi-particles 
(hadrons and/or unbound quarks and gluons) and s stands 
for their internal degrees of freedom. Here Uj(p) is the 
density-dependent mean-field acting on the quasi-particle 
j described by the field operator tpj with rrij being the 
current mass of quarks and gluons or the free mass of 
hadrons. Applying the density-dependent Hamiltonian (Tj"J 
in the partition function requires some additional con- 
straints that are needed to fulfill the thermodynamic con- 
sistency condition 29 : 



x dT' 



0, 



,dH, 
"d Pj < 







(2) 



where (A) denotes the average value of the operator A 
over the statistical ensemble. With the Hamiltonian JJJ 
the conditions (J2J can be also expressed as [3""] 



8Ui dC n v-^ dU* dC 



It can be shown 20,27,28 that the conditions (0) are 
satisfied only if the mean field Uj(p) and the correcting 
function C(p) are temperature independent. 

In the following, we consider the basic structure of the 
effective Hamiltonian (Tj"J to model the EoS of hadronic 
and quark-gluon plasma phase. 



2.1 The hadronic phase 

The hadronic phase is considered as a gas of hadrons and 
resonances in the thermodynamic equilibrium. In general, 
the particle density of species j is obtained from 

rij = rij (T, p,j -Uj) = 

d kk 2 f^k.r.,,, r ;i , (4) 

where the one-particle distribution function with an argu- 
ment z is 



fj(k,T, z) = 



exp 



'k 2 + rrij — z 



(5) 



with Cj = +1 for fcrmions and £j = —I for bosons, while 
dj is the spin-isospin degeneracy factor. The chemical po- 
tential pj is related to the baryon (pb) and strangeness 
(ps) chemical potentials 



pj = bj p B + sj Ps , 



(G) 
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with bj and Sj being the baryon number and strangeness 
of the particle j. The hadronic potential Uj = U 
described by a non-linear mean-field model |30| 



j 



is 



U j = 9r, 3 <Pl(x) + 9a J ¥>2(y) 



(7) 



where g r j > and g a j < are repulsive and attractive 
coupling constants, respectively. The effect of interactions 
results also in an additional density-dependent term C{p) 
that contributes to the thermodynamic pressure and en- 
ergy densities. If the particle interaction is taken in form of 
, the thermodynamic consistency implies that the func- 
tions <fi(x) and <p 2 (y) depend only on particle densities. 
In Ref. PHI these functions were chosen such that 



(8) 



where 



V = 



E 



9a, i 



with &i and b 2 being free parameters. The ip 2 term is in- 
troduced to get a slower than linear increase of attraction 
with density at a high compression as it happens in the rel- 
ativistic mean-field models. Having in mind that the the 
hadronic EoS will be compared with that of the quark- 
gluon plasma, it is convenient to rewrite JSJ in terms of a 
number of constituent quarks and antiquarks z/$ : 



Pj = v 3 n i = 'V'.r '-/'j ~ Uj) 



(9) 



In the original paper j.'iOj . the hadronic phase was mod- 
elled as a mixture of nucleons and A's (i.e. j — N,A). 
Following |2Jj], we generalize this approach by including 
all hadrons and resonances with the mass up to 1.6 GeV. 
One also assumes that all coupling constants scale with 
the number of constituent quarks |20| : 

= is j ([vi(P (h) T + &(P (h) )) , (10) 
where tpi and if2 satisfy Eq. Jjyl in the following form 



-C 2 (P2 - C 3 if 2 



(11) 



with pW = Ei"iPi = 3£ 



B 3 



2E 



M 3 



As compared 

to Eq. (jHJ a free parameter a is also introduced in Eq. 
(|10|> . This parameter is used to control the strength of 
the repulsive interactions at high density |27U28j . The 
parameters in Eq. Ijlljl are expressed as |27| : 



Cl 



c-2 



i9a,j/Vj) 2 ' 



C3 



b\b 2 



(9a,j/Vj) 4 



and are fixed by requiring that the properties of the ground 
state (T = and tib = n s» 0.15 fm -3 ) of the nuclear 
matter are reproduced: zero pressure, binding energy per 
nucleon of -16 MeV and incompressibility of 210 MeV. 

Solving the cubic equation (|lll) . one gets the interac- 
tion potential as 



U[ h) = v, 



Cl 



(12) 



where the function F depends on the density of quarks 
bounded inside hadrons as follows 



12 1 / 3 It / 1? \ 3 

F(t) = — rj - 2j3 V - x with v = l- + ] jf33 + -\ (13) 

Here a, (3 are proportional to the coefficients of the equa- 
tion JIIJ}: a = c 3 /9 and [3 = C2/I2 1 / 3 . 

In this representation we obtained for the hadronic 
pressure 



pt H >(T,H-Ur>) 



E d 3 



and for the energy density 



(14) 

fjik^Hj-Uf*) k 2 dk + C( P (h *>) 



where the function C is obtained from 
1 a 



(15) 



C(p {h) ) 



ci a + 1 



_ p Wj?(pW)+ (16) 



F(t) dt. 



100 



E 10 



> 

CD 




n B /n 

Fig. 1. Pressure as a function of baryon density for T = 
0. Grey and black solid lines are calculated for the modified 
Zimanyi model and ideal gas EoS, respectively. The shaded 
region corresponds to the Danielewicz et. al. constraint |21|. 



As shown in Fig. ^ the above hadronic EoS satisfies 
the constraint resulting from the nucleon flow analysis of 
heavy ion collisions in the energy range < 10 AGeV. The 
upper boundary of the shaded area is consistent with the 
constraint coming from the analysis of the neutron star 
properties |23| . In the high temperature regime there is 
also a reasonable agreement of our model with the ther- 
modynamics of the interacting pion gas from [3Tj (see 

Fig. m 
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Fig. 2. Temperature dependence of the reduced pressure and energy density for an interacting pion gas (n + p system). The 
solid line is our result, dash-dotted and dashed ones are the interacting |2J and ideal pion gas, respectively. 



2.2 The two-phase bag model 

In the MIT bag-like model, the deconfinement phase tran- 
sition is determined by matching the EoS of an ideal rel- 
ativistic gas of hadrons and resonances to that of an ideal 
gas of quarks and gluons. In the following we consider the 
two-phase (2P) model that accounts for interactions sep- 
arately in the hadronic and quark-gluon plasma phase. 
The hadronic phase is described within the phenomeno- 
logical mean-field model introduced in the previous Sec- 
tion. Following Eq. (0J, the total baryon density and the 
strangeness density in the hadronic phase can be expressed 
as 



nf = s i n i( T >Pj - u j' l) ) > 
jeh 



(17) 
(18) 



where the sum is taken over all hadrons and resonances. 
Similarly, the pressure and energy density of the species j 
are given by Eqs. Ijl4|) and (|15|) . 

In the quasi-particle approximation, the QGP phase is 
commonly described as a gas of partons (non-interacting 
point-like quarks, antiquarks and gluons) confined in a 
"bag". The non-perturbative effects associated with con- 
finement are presented by the constant vacuum energy B. 
The recent LGT results show that such an approach is 
not adequate as the EoS differs from the asymptotic ideal 
gas values even at temperatures as high as 100 T c |32| . 
The perturbative QCD results can be, however, improved 
through the so-called Hard Thermal Loop (HTL) expan- 
sion. According to the HTL perturbative expansion, the 
QCD thermodynamics at high temperature is controlled 
by quasi-particles with a temperature dependent mass 
m q (T). For /ig — one gets |32*I33] : 



2/rri\ 2 Ng 2 2 

m q (T) - m q0 = yTjj-T g 



(19) 



To model the HTL results within the mean-field ap- 
proach one introduces the quark and gluon potentials to 



reproduce the behavior of the HTL masses (|19f) in the high 
temperature limit. This in general results in an additional 
equation for the unknown gluon density. To simplify the 
problem we modify the potential so that it coincides only 
with the HTL expression for quarks. In the high temper- 
ature limit and having in mind that p ~ T 3 , the simplest 
phcnomcnological choice of the potential is 



UiPD = B (p M)i/3 i 



(20) 



where B is obtained comparing the asymptotic expansion 
of (HI with the HTP result 



B 



Ng 

16N a 



(^(2d g + 3N f d q ) 



1/3 



(21) 



with d q and d g being the degeneracy factors for quarks 
and gluons, respectively. For N c — 3 and N g = 8 one gets 



B(Nf = 3) = 0.2351 g , 
B(Nf = 2) = 0.2542 g 



(22) 
(23) 



where the strong interaction coupling constant g is 
treated as a free parameter. 

The thermodynamic self-consistency conditions re- 
quire that the mean-field contribution to the pressure and 
energy density in equations like Eqs. i|14fl and (|15fl is re- 
spectively 



(24) 
(25) 



where the plasma particle density p^ p ' = X^'eg q g Pj an d 
the bag constant B is included in the correcting function 
C. 

With such mean-field potentials the pressure and en- 
ergy density in the plasma phase carried by u, d and s 
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quarks and antiquarks is obtained as 



q 

-C(p (pl) ) , 



(26) 



(27) 



To quantify these observables we use the quark masses 
m u = md = 65 MeV and m s = 135 MeV, the gluon mass 
nig = 700 MeV and the bag constant B 1 / 4 = 207 MeV. 
Such parameters yield a transition temperature T c ss 
170 MeV in agreement with the recent lattice result ob- 
tained for the vanishing net baryon number |2] . For mass- 
less gluons the equation of state has a simple form 

P 9 (T) = ^T\ £g (T) = 3p g (T) = i£ T 4 (28) 



90 



30 



with d g =16 



The baryon number and strangeness density in the 
quark-gluon plasma are obtained following Eqs. (|17|) and 
iJTSJl from 



£ bj rij(T,fij — Uj) , 



= £ Sj rij(T,fij — Uj) 



(29) 
(30) 



j£g,q,q 



The equilibrium between the plasma and the hadronic 
phase is determined by the Gibbs conditions for the ther- 
mal (T^ = T H ), mechanical (jfi — p H ) and chemical 
(fig = fig, fig = fig) equilibrium. At a given temper- 
ature T and baryon chemical potential fiB the strange 
chemical potential fis is obtained by requiring that the net 
strangeness of the total system vanishes. Consequently, 
the phase equilibrium condition and strangeness conser- 
vation imply that: 



p H (T,p j -U^)=pQ(T,^-U^) 



n B 



(l-A)n£(T,^--[/ (pi) ) + 



+\n Q B (T,p j -U { ? l \ 



0=(l-X)n§(T, H -U^) + 



(31) 
(32) 

(33) 



+Xn^(T,p j -U^), 



where A = Vq/V is the fraction of the volume occupied by 
the plasma phase. The phase boundaries of the coexistence 
region are found by putting A = for the hadron phase 
boundary and A = 1 for the plasma boundary. By con- 
struction the 2P EoS results in the first-order phase tran- 
sition with discontinuous behavior of energy and baryon 
densities. 

According to the Gibbs condition [2], the number 
of thermodynamic degrees of freedom that may be var- 
ied without destroying the equilibrium of a mixture of 



r phases with n c conserved charges is Af — n c + 2 — r. 
For the considered hadron-quark deconfincmcnt transi- 
tion r = 2. If the baryon number is the only conserved 
quantity then n c = 1 and Af = 1. Thus, the phase bound- 
ary is one-dimensional, i.e. a line. The Maxwell construc- 
tion for the first-order phase transition corresponds to 
r = 2 and n c — 1. When both the baryon number and 
strangeness are conserved, that is when n c = 2, one has 
Af = 2 and therefore the phase boundary is a surface. 
In such a system, a standard Maxwell construction is no 
longer possible pSllSrJ] 1 . 

When two phases coexist, the system is in general not 
homogeneous because the phases occupy separate domains 
in space. We do not, however, explicitly account for such 
a domain structure or a possible surface energy contri- 
bution to the equation of state. The only consequence of 
the phase separation in the considered 2P model is that 
the interaction between quasi-particles in the plasma and 
hadronic phase are neglected. This is different from the 
statistical mixed phase model that will be discussed in 
the next subsection. 

The resulting phase boundaries between the hadronic 
phase and the quark-gluon plasma in the 2P model are 
shown in Fig. At T — the coexisting region appears at 
n-B /no ~ 8. This density is by factor two larger than that 
one obtained in the conventional MIT-bag like model (20] 
(see also Ref. [22] )• This is because in our calculations the 
quarks and gluons are treated as massive quasi-particles. 
As it will be shown in the next Section, the finite mass of 
quasi-particles is needed in the quark gluon plasma to get 
the EoS that is consistent with LGT results. 



2.3 The mixed-phase model 

In the 2P model the interactions between quark, gluons 
and hadrons are entirely neglected in the coexistence re- 
gion. In the following we introduce the MP model where 
such interactions are possible. The underlying assumption 
of the MP model [2*7ll2*H] is that unbound quarks and glu- 
ons may coexist with hadrons forming a homogeneous 
quark/gluon-hadron phase. Since the mean distance be- 
tween hadrons and quarks or gluons in this mixed phase 
may be of the same order as that between hadrons, the in- 
teraction between all these constituents (unbound quarks, 
gluons and hadrons) plays an important role as it defines 
the order of the phase transition. 

Under a quite general requirement for the confinement 
of color charges, the mean-field potential of quarks and 
gluons in the plasma phase is approximated as 



U q (pW) = U g (pW) 



(p'h _ 



A 



(p( P i))i 



(34) 
7>0 , 



1 In |35I36) the baryonic and electric charge conservation was 
considered in application to a nuclear liquid-gas phase transi- 
tion. As to strangeness conservation the emphasis was made 
mainly on the strangeness distillation effect 37 . Phase bound- 
aries for this case were studied in detail in |38| . More complete 
list of appropriate references can be found in HI]. 
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Fig. 3. The phase boundary calculated in the 2P model with the physical values of parameters as explained in the next Section. 



500 1000 1500 2000 2500 
MeV 



where p 



(pi) 



The second term in Eq. (|35fl is 



introduced to account for the growth of the quasi-particle 
mass with the density as that obtained in the HTL ap- 
proximation (see Eq. <\2A\). The first term in Eq. (|35() re- 
flects two important limits of the QCD interactions. For 
p(pl) _> o, this potential term approaches the infinity, i.e. 
an infinite energy is necessary to create an isolated quark 
or gluon that corresponds to the confinement of color ob- 
jects. The other extreme limit of infinite density is consis- 
tent with the asymptotic freedom. 

The generalization of the mean-field potential from 
Eq. (|35[l to the case of the mixed quark-hadron phase is 
obtained replacing in Eq. I|35|) by the total density of 
quarks and gluons p( mp ) with 



P 



(mp) 



Pq+ Pq+ Pg+VY1 



iPj 



,(PQ 



VP 



(/,.) 



(35) 



The presence of the total number density p( mp > in Eq. i|35fl 
implies interactions between all components of the mixed 
phase. For r\ = there is no interaction between hadrons 
and unbound quarks and gluons. This case corresponds 
to such a strong binding of hadron constituents that the 
presence of free color charges in their surrounding does not 
result in their color polarization, i.e. hadrons remain color 
neutral and do not see quarks and gluons outside hadron. 
Thermodynamically, the potential with r\ = implies the 
first order phase transition. For r\ = 1 there is a very 
strong color polarization of hadrons. Consequently, there 
is no difference between bound and unbound quarks and 
gluons. This approximation was used in |27U28) . Here we 
consider 77 as a free parameter that is chosen in a way to 
reproduce the LGT results for the QCD equation of state. 

The hadronic potential in the Hamiltonian Q was de- 
scribed by a non-linear mean-field model. However, the 
presence of unbound quarks and gluons will modify this 
hadronic interaction due to the polarization of color 
charges. Thus, in general 



U 



(mp) 



u 



(h) 



u 



(h-pl) 



J 



(36) 



The constraints imposed by the thermodynamic con- 
sistency conditions © can be used to find the potential for 



the interaction of unbound quarks/gluons with hadrons as 

EM 



u. 



(h-pl) 



A 



A 



(p("ip))i (r]p( h ))~f 

B [(p m P)) 1 /3„( 7?/ ,W)l/3- 



(37) 



Consequently, the pressure and the energy density in 
the MP model are obtained from 



MP 



Y^PdT,P,-U { 3 mp) )-C(p^), 
jeh 



P 



(38) 



(mp) 



) + 



(39) 



where 



C{p( mp) ) 



p (.h) 
jA 

B 

4 



_ (p w r+ i 

F{ P {h) ) - 



(P 





(roph 1-7 



F{t)dt - 



(pW ) 4 / 3 - ( W (,i) ) 4/E 



(40) 



The MP model described above exhibits a crossover 
deconfmcment phase transition. The transition tempera- 
ture T c corresponds to a maximum in the T-dependence 
of the heat capacity at the given value of [Lb (see the next 
Section). The resulting phase boundary is shown in Fig. 
0] At T < 50 MeV the maximum of the heat capacity is 
not well denned. The calculation in Fig. 0] was performed 
with the physical values of the parameters as introduced 
in the next Section. 
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Table 1. Parameters of the interpolation formulae 1(41(1 
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Fig. 4. The phase boundary calculated in the MP model 
(solid line) . The dotted and dot-dashed lines correspond to the 
state where with the fraction of unbound quarks is 0.5 and 0.6, 
respectively. 



In the MP model hadrons survive at T > T c . If the 
fraction of unbound quarks is defined as p^ / (p( pl > +p^), 
then one can see from Fig. 0] that at ps = and the 
temperature as high as T ~ T c ~ 100 MeV there is still 
40% quarks that are bounded inside hadrons. 



3 The comparison with the lattice data 

The Lattice Gauge theory is the only approach that allows 
to extract the physical EoS of QCD medium. To further 
constraint the phcnomcnological models for the EoS intro- 
duced in the last Section, we will compare their predictions 
with the available LGT results. We focus mainly on the 
recent LGT findings obtained in (2+l)-flavor QCD at the 
finite temperature and chemical potential 

In order to use the mixed phase and 2P models for the 
further comparison with lattice results one needs, however, 
to take into account that lattice calculations are generally 
performed with quark masses heavier than those realized 
in the nature. Consequently, the hadron mass spectrum 
generated on the lattice is modified. 

In Rcfs. 4,5 the ratio of the pion mass to the mass 
of the p meson is around 0.5-0.75, which is roughly 3 times 
larger than its physical value. Thus, to compare the model 
predictions with LGT results the hadron mass spectrum 
used in the model calculations should be properly scaled. 
For this we use a phenomenological parametrization of the 
quark mass dependence of the hadron masses rrij (x) that 
was shown in Refs. |15L 16 : to be consistent with the MIT 
bag model results as well as with LGT findings. For the 
non-strange hadrons this parametrization reads |15Ul(j| : 



mj{x) 



VljCLiX 



m. 



1 + d2X + CI3X 2 + d4,X 3 



a>5X" 



.(41) 



Here x = m^/y/a, V[j is the number of light quarks inside 
the non-strange hadron ( i.e. z/y = 2 for mesons and vij = 
3 for baryons) and a = (0.42 GeVf . 



Oi 


a-2 


a-3 


0,4 15 


0.51 




0.115 


-0.0223 0.0028 




rrij 







For strange hadrons that carries strangeness Sj = 1 
and Sj — 2 we have, respectively 



m(x) 1.7-0.42 3i 



m,j{x) 



0.5788a; - 



(1 + 0.068x) ' 
0.42^= 
(0.4758 + 0.0142a;) ' 



(42) 
(43) 



Simultaneously with the change of the hadron mass spec- 
trum with the pion mass one needs to account for the 
shift of the transition temperature T c with m^. We use 
the parametrization that is extracted from LGT calcula- 
tions |T5). 



T, 



0.4 + 0.04(1) ^ 
' la 



(44) 



To compare our phenomenological model EoS with 
that obtained on the lattice in Refs. we use the mod- 
ified hadron mass spectrum form Eqs. 1)4 l|l - 143 (1 corre- 
sponding to the pion ~ 508 MeV as fixed in these LGT 
calculations. In the deconfincd phase the current quark 
masses and gluon mass are in general also free parame- 
ters. In the present calculations we fixed m u — rrid = 65 
MeV, m s = 2.08 m u and m g ~ 700 MeV as followed from 
the successful description of the quark sector of the above 
LGT data in terms of the quasi-particle model [8]. 

We look for a phase transition at the appropriate tem- 
perature T c defined by Eq. I|44|l by varying mainly the bag 
constant B in the 2P model or strength parameter A in 
the mixed phase model. The further fine tuning is carried 
out by means of remaining parameters (one parameter in 
the 2P model and three ones for the mixed phase model) 
to get the best description of LGT findings on tempera- 
ture dependence of different thermodynamical quantities. 

In the 2P model, where the two phases do not interact 
with each other, the critical temperature T c is governed 
mainly by the value of the bag constant B and the parame- 
ter a that characterizes the hardness of the EoS. Choosing 
B X I A = 223 MeV and a = 2.1 one gets T c w 176 MeV and 
e/T 4 \ Tc = 7.84 to be consistent with the LGT results. We 
have to stress, however, that for some values of the param- 
eters, e.g. for too heavy masses, the set of Eqs. 1(31(1 - ((3411 
may have no solution. 

In the MP model the critical temperature is defined at 
the position of the maximum of the heat capacity 

c v = de/dT\ v=const 

The value of T c depends mainly on parameters that quan- 
tify the quark/gluon interactions. With yl 1 /( 3 T+ 1 ) = 270 
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MeV and 7 = 0.3 the accepted value of the critical tem- 
perature is seen in Fig. El to be 188 MeV. As was noted 
in Ref. j2U], 7 = 1/3 corresponds to a string-like quark 
interaction. 



120 r- 
100 - 




T, MeV 

Fig. 5. Temperature dependence of the reduced heat capacity 
at hb = 

In Figs. El and we show the comparison of the MP 
and the 2P model predictions with LGT data obtained for 
the thermodynamic pressure p/T 4 and the energy density 
e/T 4 at the finite T but for fi B = 0. 

The lattice calculations in Refs. 0]E] were done on the 
lattices with N t — 4 temporal extension. To account for 
the finite size effects, the LGT results have to be extrap- 
olated to the continuum limit corresponding to Nt — > 00. 
In general, such a procedure requires a detailed LGT cal- 
culations on the lattices with different N T . In Refs. 0]E] to 
account approximately for the finite size effect the N T = 4 
data for the basic thermodynamic quantities were cor- 
rected being multiplied by the constant factors: cq — 0.518 
and = 0.446 for fi B = and ps 7^ respectively. These 
factors were determined from the ratios of the Stefan- 
Boltzmann ideal-gas limit for the thermodynamic pressure 
to its corresponding values calculated on the lattice with 
N t = 4. 

As it is seen in Fig. El a smooth T-dependence of 
pressure in a deconfined phase may be quite well repro- 
duced within both the MP and 2P model. However, in the 
hadronic phase, that is for T/T c < 1, the models overes- 
timate LGT results from Refs. |4"|y>]. In Fig. also shown 
are LGT results for (2 + l)-flavor QCD at fis — from 
the Bielefeld group jSH]- Improved gauge and staggered 
fermion actions were used there on the lattices with tem- 
poral extent of Nt = 4 and Nt = 6. These data were 
also extrapolated to the chiral limit [HHI- As seen in Fig. 
E| the Bielefeld data exhibit a smaller limiting pressure 
as compared to 0]E] and essentially higher pressure in 
the hadronic sector, though the pion mass is = 770 
MeV in the latter calculations. Our models are seen in 
Fig. El to coincide with Bielefeld results in the confined 
phase. The strongly suppressed pressure at T < T c found 
in Refs. |3]E] is non-physical and could be partly related 
with too simplified procedure to extrapolate LGT results 



to a continuum limit when applying the same constant 
scaling factor for all values of temperatures. 

The energy density shown in Fig. behaves differ- 
ently in the MP and 2P model. As it is expected in the 
2P model, that exhibits the first order phase transition, 
the e/T 4 suffers a jump at the critical temperature. This 
jump corresponds to the energy density change by As ~ 
0.9 GeV/fm 3 . The LGT results on the temperature depen- 
dence of e/T 4 are seen in Fig. to be noticeably better 
reproduced within the MP than with the 2P model. This is 
because the MP model exhibits a crossover type transition 
as also found in the above LGT calculations. The differ- 
ence between LGT results obtained with an improved and 
standard action is also seen on the level of energy density. 

Having established the model parameters at ps = 
we can further study the model comparisons with LGT 
results at the finite baryon density. The temperature de- 
pendence of pressure and energy density for finite values 
of [Ib is shown in Figs.ElandElin terms of the "net bryonic 
pressure" Ap/T 4 = (p(T, p B ) ~ p(T, p B = 0))/T 4 and the 
"interaction measure" A/T 4 = (e - 3p)/T 4 . 

The T dependence of Ap/T 4 for different values of [is 
is quite well reproduced by both the MP and 2P model. 
The fall of Ap/T 4 for T > T c is entirely determined by 
the value of the coupling g that describes the strength of 
the interactions of quasi-particles and their effective mass. 
The observed fall does not require any artificial reduction 
of the number of quark-gluon degrees of freedom. It turns 
out that in both models a similar value of g — 0.5 is 
necessary to reproduce LGT results. 

The interaction measure A/T 4 exhibits a rather sharp 
maximum slightly above T c with the shape of In- 
dependence that is weakly changing with ps- In general, 
both models reproduce the above properties of the inter- 
action measure. However, quantitatively the A/T 4 is over- 
estimated in the 2P model and underestimated in the MP 
model near the maximum. 

The interaction measure characterizes the strength of 
interactions in a system. It is equal to zero for the EoS of 
the ultrarelativistic ideal gas of massless particles where 
e = 3p. In the considered 2P model and at T > T c we are 
dealing with a gas of massive quarks and gluons interact- 
ing via the HTL-like potential l)24|l. In contrast, in the MP 
model at T > T c , there are interacting unbound quarks, 
gluons and bound quarks within hadrons. The fraction of 
the bound quarks amounts in about 85 % at T ~ T c (which 
allows to describe this region in terms of the resonance gas 
model |15Ulfi| ) and almost vanishes at 3T C (~ 5%). In this 
context the quark-gluon plasma may be considered as a 
strongly interacting correlated system 013 ■ ^ n the con- 
fined phase there is an admixture of quarks at T < T c 
until about 0.9 T c . This property of the model is very es- 
sential for a possible explanation of the "horn" structure 
in the K + /ir + excitation function due to manifesta- 
tion of the strangeness distillation effect near the critical 
end point |4~2*] . 

The model comparison with LGT results for the 
baryon density is shown in Fig. 1101 It is clear from this fig- 
ure that in the hadronic sector the baryon density ns/T 3 
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Fig. 6. The reduced pressure at \ib = in 2P (the left panel) and MP (the right panel) models. Circles are the lattice data 
for the (2+l)-flavor QCD system |4I5) multiplied by co, squires are the Bielefeld group data for the same case |39) (see also 
results cited in |15|.^ 





■ ■ ■ ■ ■ i ■ ■ ■ ■ i ■ ■ ■ ■ i 



0.5 1.0 1.5 2.0 2.5 3.0 

T/T c 

Fig. 7. The reduced energy density at /is = in 2P (left panel) and MP (right panel) models. Notation is the same as in Fig. 

m 



obtained on the lattice is smaller than that predicted by 
the models. However, above T c there is quite a good agree- 
ment of the model with LGT results. This is particularly 
a case for the MP model which shows a better description 
of LGT data near the phase transition. 

In our models the absolute values of Ap/T 4 , A/T 4 and 
ns/T 3 are strongly affected by the parameter 77 appearing 
in Eq. Ij35|l . In the actual calculations the r\ = 0.025, thus 
it is essentially smaller than r\ — 1 that was found in 
our earlier parametrization based only on the LGT data 
obtained for [Ib = 20 . If 77 — 1 is to be substituted 
in our actual calculations, then all the above quantities 
would increase by a factor of two. 

The properties and the behavior of the LGT thermo- 
dynamics at finite T and \ib have been recently discussed 
in the context of the Polyakov-loop-extended Nambu and 
Jona-Lasinio (PNJL) model ^21- This PNJL model repre- 
sents a minimal synthesis of the spontaneous chiral sym- 
metry breaking and confinement. The model correctly de- 
scribes the pion properties but obviously is not applicable 
near the nuclear ground state. It also does not contain the 



resonance contributions to the QCD thermodynamics nor 
the hadronic correlations below and above T c that are es- 
sential near the phase transition. Nevertheless, the PNJL 
model reproduces the LGT data [H] obtained in 2 flavor 
QCD on the pressure difference and the quark number 
density at various temperatures and chemical potentials 
remarkably well ^3|- However, in the PNJL model the in- 
teraction measure Z\/T 4 was found to be underestimated 
by ~ 25% similarly as seen in FigUHfrom our MP model 
comparison with the (2+l)-flavor QCD results obtained 
in LGT. 

The phenomenological models considered here describe 
the EoS of the QCD matter in a broad range of thermal 
parameters that includes the hadronic and quark gluon 
plasma phase. These models are also applicable in cold 
nuclear matter as they satisfy essential phenomenological 
constraints expected near nuclear saturation. In heavy ion 
collisions, dense QCD matter created in the initial stage is 
expected to thermalize and expand without further gener- 
ation of the entropy S. In the realistic expansion scenario 
some particles may be crated and/or absorbed implying 
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T/T c 



Fig. 8. Temperature dependence of the reduced pressure 

410 and 530 MeV (from the bottom) within 2P ( the left pane! 

the (2+l)-flavor system 14151 multiplied by c M . 



1.0 



0.8 - 




T/T c 

) — p(0))/T 4 at the baryon chemical potential \ib = 210, 330, 
and MP (the right panel) models. Points are lattice data for 




changes in the total entropy of the system. In general, it is 
more convenient to consider the EoS at fixed entropy per 
baryon (S/Nb)- This thermodynamic quantity should be 
strictly conserved in an equilibrium case and is also less 
affected by any possible particle losses or creation during 
the expansion stage. The predictions of our models for the 
evolution path in the (T, /Lts)-plane as obtained from the 
condition of fixed S/Nb is shown in Fig. EI There are still 
no such isentropic lattice data for (2+l)-flavor system. 
Recently, the isentropic EoS was obtained on the lattice 
for 2-flavor QCD at finite ^b |43| . however still for non- 
physical mass spectrum that corresponds to the pion mass 
~ 770 MeV. These data are plotted in Fig.lllltogether 
with our model results obtained with the EoS parameters 
that are fixed for ~ 508 MeV and for (2+l)-flavor 
system. 

In general, in the high-temperature deconfined phase, 
one should not expect a large difference between 2 and 
(2+l)-flavor thermodynamics. The value of the quark 
mass in the quark-gluon plasma is also not relevant ther- 



modynamically if m q /T < 1. In the hadronic phase the 
number of quark flavors is as well not essential and leads 
only to a moderate change of the global thermodynamics. 
However, here the value of the quark/pion mass is of par- 
ticular importance as it influences the hadronic mass spec- 
trum. Due to the non-physical and still large pion mass 
used in the actual LGT calculations it is not straightfor- 
ward to associate the values of the reduced entropy with 
the specific bombarding energy. In particular, as noted 
in Ref. 03| j the correspondence of S/Nb — 30, 45 and 
300 to the AGS, SPS and RHIC energies, respectively, is 
only a rough approximation. The QGSM transport model 
results 03] for central Pb+Pb collisions at the top SPS en- 
ergy show that the isentropic regime is reached after about 
1 fm/c with S/Nb ~ 25. Also calculations performed in 
terms of 3-fluid relativistic hydrodynamic model show 
that the isentropic expansion of central Pb+Pb collisions 
at the bombarding energy 158 and 30 AGeV results in 
S/Nb ~ 30 and 15 0B] respectively. Thus, the above dy- 
namical models imply noticeably lower values of S/Nb 
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Hb> MeV 

Fig. 11. Phase trajectories in the T — \ib representation. 
Circles, squares and triangles are the lattice 2-flavor QCD 
results [H| for S/N B = 30, 45 and 300, respectively. The 
(2+l)-flavor model predictions are plotted by solid (2P EoS), 
dashed (MP EoS) and dot-dashed (Hadronic EoS) lines for ev- 
ery value of the reduced entropy. The dotted line parameterizes 
the freeze out curve I44| . 



than that obtained within actual LGT calculations |43j . 
The main origin of the above differences is related with 
still too large quark mass used on the lattice. 

As seen in Fig. El the MP and 2P models reproduce 
the general trend of the lattice trajectories. The lattice 
evolution paths are just between 2P and MP model pre- 
dictions. With increasing S/Nb these differences are no- 
ticeable smaller. The hadronic EoS predicts higher initial 
temperatures, however all three phenomcnological models 
give similar results for the freeze-out temperature. It is in- 
teresting to see in Fig. llll that irregularity appearing near 
the turning point of the lattice trajectory correlates with 
the flattening of the T-dependence in the Gibbs mixed 
phase resulting in the 2P model. 



The phenomenological model results, discussed so far, 
were obtained assuming the hadronic mass spectrum that 
corresponds to the pion mass = 508 MeV. The extrap- 
olation of the EoS to the physical limit is quite straight- 
forward. It amounts in replacing the mj{m 7r ) masses by 
their physical values. The quark and gluon masses are 
kept to be the same as being extracted from the LGT 
data. This approximation is justified since the change of 
m q in the interval 5 < m q < 70 MeV does not influ- 
ence the thermodynamics in the plasma phase |H] very 
much. Clearly, taking the physical limit in the EoS also 
requires to account for the shift in T c . In the 2P model 
the critical temperature is recalculated according to Eq. 
(14411 and fitted in the model by the bag B constant and 
the coupling g to satisfy also the condition that the crit- 
ical energy density e c /T 4 ~ 6 ± 2 as found in LGT 
|15| . Within the 2P model the physical limit is achieved 
choosing: B 1 / 4 = 207 MeV and g = 0.7 which results in 
T c = 173.3 MeV and e c /T 4 = 7.83. 

The extrapolation of the MP model EoS to the physical 
limit is less transparent due to a rather strong nonlinear 
relation between the hadronic and plasma phase. In this 
model the physical limit is approximately accounted for by 
replacing the LGT mass spectrum by its physical form. All 
further parameters that are require in the MP model to 
quantify the EoS are kept the same as that found in the 
comparison of the model predictions with the LGT results. 
With the above chosen model parameters the crossover 
deconfmement transition appears at T c = 183 MeV. Note 
that the phase boundaries in the preceding Section were 
calculated for these physical parameters of the EoS. 

4 Summary 

We have formulated two different phenomenological mod- 
els for the equation of state within the quasi-particle ap- 
proximation of the QCD matter: the two phase (2P) 
model with the first order deconfmement transition and 
the mixed phase (MP) model in which transition from 
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hadronic phase to quark-gluon plasma is of the crossover 
type. 

In our approach both the hadronic and the quark- 
gluon plasma phase are considered to be the non-ideal sys- 
tems. The interactions between constituents are included 
within the mean-field approximation. The modified mean- 
field Zimanyi model is applied to describe the interacting 
resonance gas component. In this approach, the saturation 
properties of a symmetric nuclear matter in the ground 
state are reproduced correctly and the Danielewicz con- 
straints resulting from heavy-ion collision data at inter- 
mediate energies are well fulfilled. 

The quark-gluon phase in the 2P model is con- 
structed as a massive quasi-particle system supplemented 
by the density-dependent potential term which simulates 
the HTL interactions. The first order phase transition 
from the hadronic phase to the deconfined quark-gluon 
plasma is constructed within the 2P model by means of 
the Gibbs phase equilibrium conditions. 

In the MP model the coexistence and correlations be- 
tween quarks/qluons and hadrons are assumed near de- 
confinement. In addition to the HTL-likc interaction term 
a string-like interaction is introduced between both un- 
bound quarks/qluons and quarks that are confined within 
hadrons. In this model we are dealing with strongly in- 
teracting QCD matter which exhibits a crossover-type de- 
confinement phase transition. 

The models are constructed in the way to be thermo- 
dynamically consistent and to reproduce the properties of 
the EoS as calculated on the lattice. The limited set of 
model parameters is defined from the constraints imposed 
by the recent lattice data on the temperature and chemical 
potential dependence of the basic thermodynamical ob- 
servables. The comparison of the model predictions with 
LGT data was performed within the same set of approxi- 
mations as used on the lattice. Of particular importance is 
a correct treatment of the hadronic mass spectrum which 
in the LGT calculations is non-physical due to the still 
too large value of the quark mass. 

Keeping in mind a principal difference between the 
first order and crossover type phase transition both the 
2P and MP model were shown to provide quite satisfac- 
tory description of the LGT thermodynamics for (2+1)- 
flavor QCD. Both models reproduce the T and /is depen- 
dence of the main thermodynamic quantities in a broad 
range of thermal parameters. The observed deviations of 
the model predictions from the lattice results near T c and 
in the hadronic sector for (2+l)-flavor case may be, to a 
large extend, attributed to uncertainties in the LGT data 
due to the finite size effect. The predicted isentropic tra- 
jectories in the phase diagram were shown to be consistent 
with that recently calculated on the lattice within the 2- 
flavor QCD. 

The phenomenological equations of state constructed 
here satisfy all physically relevant constraints expected in 
the cold and excited nuclear matter. These EoS can be 
applied in a broad parameter range that covers the region 
of dcconfmement transition in QCD. Thus, both the MP 
and 2P EoS could be used as an input in dynamical mod- 



els that describe the space-time dynamics and evolution 
of a medium created in heavy ion collisions. Within hy- 
drodynamic models our EoS can be important to study 
the role and influence of deconfinement and the order of 
the phase transition on physical observables. Such studies 
are in progress. 
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